Flat branches and pressure amorphization 
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After summarizing the phenomenology of pressure amorphization (PA), we present a theory of 
PA based on the notion that one or more branches of the phonon spectrum soften and flatten with 
increasing pressure. The theory expresses the anharmonic dynamics of the flat branches in terms 
of local modes, represented by lattice Wannier functions, which are in turn used to construct an 
effective Hamiltonian. When the low-pressure structure becomes metastable with respect to the 
high-pressure equilibrium phase and the relevant branches are sufficiently flat, transformation into 
an amorphous phase is shown to be kinetically favored because of the exponentially large number of 
both amorphous phases and reaction pathways. In effect, the critical-size nucleus for the first-order 
phase transition is found to be reduced to a single unit cell, or nearly so. Random nucleation into 
symmetrically equivalent local configurations characteristic of the high-pressure structure is then 
shown to overwhelm any possible domain growth, and an "amorphous" structure results. 



I. INTRODUCTION 

Since its first observation in Gd2(M2CU) in 1972,0 and 
especially in ice nearly twenty years agoJ3 amorphization 
induced by pressure (or "pressure amorphization" , PA) 
has been discovered in materials of essentially all binding 
types, from ionic to metallic, which have more than one 
atom per primitive cell. The disappearance of the X-ray 
diffraction pattern of the low-pressure crystalline phase, 
or a concomitant broadening of peaks in the vibrational 
spectra, is taken to be evidence of an emerging amor- 
phous state. These signatures typically appear within a 
pressure window whose average value and overall range 
are strongly material dependent. 

The phenomenology and simulations of materials dis- 
playing OTpSfiHFp amorphization have been thoroughly 
reviewed. eIqEHju ft is not sufficiently generally recognized 
within that extensive literature that there are two dis- 
tinct classes of materials which could be labelled as pres- 
sure amorphized. In the first, a conventional glass, the 
atomic structure is not topologically equivalent to any 
crystalline structure. In the second, random displace- 
ments of the nuclei occur which do not destroy crystalline 
topology. 

Despite extensive investigation, and in part because 
the above distinction has not been clearly made, a general 
physical picture of the microscopic processes underlying 
PA has yet to emerge. On the theoretical side, the main 
efforts have focused on two materiai^pdiibitfc 

(a-quartz)llay>BE 



ond class of PI 

and a-AlP04Jl3'E3E3£3 Both interatomic potentials and 
fully first-principles approaches have been used to study 
these systems. While it remains uncertain whether PA 
has been actually observed in a simulation, it seems to 
be generally accepted that there is a correlation between 
pressure amorphization and the softening (i.e., phonon 
frequencies approaching zero) of a nearly dispersionless, 
or flat, acoustic phonon branch. In addition, it has 
been suggested that the occurrence of an amorphous 
phase sbojiki be strongly related to the dynamics of the 
system-BO 



In this paper, we discuss an explicit model for mate- 
rials exhibiting the second class of PA. On the basis of 
inferences drawn from both phenomenology and simu- 
lations, we first argue that these materials share many 
common features and then show how to construct an ef- 
fective Hamiltonian that incorporates these features. As 
we demonstrate, this microscopic model is simple enough 
to analyze in detail yet still relevant to real materials. 

The paper is organized as follows. In Section II we 
discuss the elements of our theory of PA on the basis of 
both phenomenology and simulations. In Section III wc 
review some elements of lattice vibration theory. In Sec- 
tion IV, we present an effective Hamiltonian suitable for 
the discussion of the physics underlying PA, and reduce it 
to the simplest model capable of incorporating the essen- 
tial features discussed in Section II. In Section V, we de- 
scribe how PA is manifested in the simple model through 
an analysis of the topology of the model's potential en- 
ergy surfaces. Sections VI, VII, and VIII address the 
conditions for the occurrence of pressure amorphization 
in the 2-site, 3-site, and iV-site models respectively. In 
Section IX, we discuss the relation of our considerations 
to conventional nucleation and growth models of struc- 
tural transitions. We summarize our results and draw 
conclusions in Section X. 



II. PHENOMENOLOGY AND SIMULATIONS 



Many transition patterns,,-have been observed dur- 
ing pressure amorphizationfltM Of these, we select five 
(types I-V) for explicit discussion: 

Type I materials transform from a crystalline phase C 
to an amorphous phase a at pressure p — p ca . Then, 
after p is reduced from above p ca , they remain amor- 
phous down to a lower pressure p aci below which they 
revert to C, as is the case, for example, in AIPO4.E3 (The 
hysteresis p ca -p a c is large, typically of the order of p ca it- 
self.) Materials of type II remain amorphous upon release 
of pressure, making possible detailed examination of the 
properties of the amorphous phase. This appears to be 
the case for amorphous a-Si02 , which has been found to 



2 



retain elastic anisotropyBdO Another example is crys- 
talline Fe 2 Si04-fayalitc, which is antiferromagnetic; the 
pressure- amorphized fayalite remains antiferromagni 
with the Neel temperature essentially unchanged 
from which we infer that the a- and C-phase topologies 
are the same. Type III materials become amorphous 
above p ca , later transforming to a new crystal structure 
C when p is held steady above p ca for sufficiently long. 
Materials of type IV behave the same as those of type I as 
p is reduced from above p ca , but they transform to a new 
crystal structure C from the amorphous phase a when 
pressure is increased beyond p ac i > p ca . Finally, materi- 
als of type V transform to C at p QC < upon decompression 
from above p ac , finally reverting to C below p cc i. 

As already mentioned in the Introduction, lattice dy- 
namical calculations have shown a considerable softening 
of the entire lowest band of a-quartz, with a soft mode 
instability occurring first at the zone boundary at pres- 
sures at which PA is observedEiEj or at somewhat higher 
pressuresJla Also, it has been argued that the transition 
mechanism in a-AlPC>4 is similar to that of a-quartz. Ej'l£| 

From the phenomenology and calculations described 
above, the following picture of PA emerges: The ran- 
dom displacements in the amorphous phase are associ- 
ated with the incipient instability of at least one nearly 
flat phonon branch. The amorphous phase is the result 
of a kinetically hindered transition between two closely- 
related crytalline phases C and C , where C is a low- 
pressure crystalline phase and C is a high-pressure one. 
The transition from C to C is expected to be strongly 
first order, and also displacive, meaning that the atomic 
structure of a remains topologically equivalent to that 
of C. (This latter feature is to be contrasted with 
the significant changes in topology that usually accom- 
pany ordinary glass formation.) In the amorphous phase 
a, random displacements are large enough to eliminate 
the diffraction pattern through reduction of the Debye- 
Waller factor, but not so large so as to disrupt the topol- 
ogy of the crystal structure associated with C. 



III. NEARLY FLAT BRANCHES; LOCAL 
MODES AS LATTICE WANNIER FUNCTIONS 

Let us assume that we can work in the Born- 
Oppenheimer approximation, supressing all the elec- 
tronic degress of freedom. The Hamiltonian is, then, 



H[P,R] =T[P] + V[R], 



(1) 



where P and R are the sets of nuclear momenta and 
positions. T[P] is the nuclear kinetic energy, 



p2 

TIP] = V— , 



(2) 



where Pi is the momentum of the z-th nucleus and Mj is 
its mass. V[R] is the nuclear potential energy 



V[R] = E[R] + V NN [R], 



(3) 



where E[R] is the electronic ground state energy and 
Vnn[R] the internuclear interaction energy for configu- 
ration R. 

Let Rq be a set of positions corresponding to a local 
minimum of the energy that we choose as a reference (e.g. 
the low-pressure phase C). Suppose that R differs from 
Ro through the uniform strain tensor e and the nuclear 
displacements u'\ 



R = R Q -(l + e) + u', 



(4) 



where all site, vector, and tensor indices are understood, 
and the dot product implies a sum over all relevant site 
and component indices. Expanding V to second degree 
in e and u', we get: 



V[R] = V[R \ + -u 



1 , d 2 V[R ] 



, d 2 V[R ] 1 d 2 V[R Q ] 



Making the substitution v! = u + uq with 

-l 



u 



d 2 V 
du'du' 



d 2 V 
du'de 



we obtain 



H 



T + -u- 



d 2 V 



where 



C = 



d 2 V d 2 V ( d 2 V 



dede dedu' \ du'du' 



d 2 v 

du'de 



(6) 



(7) 



(8) 



Thus, when d 2 V/ du'du' becomes soft, that is, when it 
has eigenvalues which are nearly zero, C will also be- 
come soft according to Eq.(||), provided the correspond- 
ing eigenfunctions of d 2 V /du'du' couple to e. In mate- 
rials exhibiting PA, the modes which soften under pres- 
sure will couple to strains associated, with elastic moduli 
that likewise soften under pressure.0 We can thus focus 
primarily on the softening of d 2 V/ du'du' and treat the 
strains implicitly, recognizing that the u in Eq.(7) con- 
tains them. 

The transformation 



(9) 



la 



diagonalizes H. In Eq.(||), the site index of Eq.(||) has 
been decomposed into a unit cell index I, with Ri the 
center of the Z-th cell, and an internal index a. A is the 
branch index of the normal mode fcA, k its wave vector in 
the first Brillouin zone, and S\ a (k) is a unitary matrix 
for each fc. The squares of the corresponding normal- 
mode frequencies u>k\ are eigenvalues of the dynamical 
matrix M' 1 ' 2 ■ d 2 V/du'du' ■ M" 1 / 2 . Thus if uj kx has 
a soft, flat branch in which ujkx — * for much of the 
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Brillouin zone, so does d 2 V / du' du' , as the mass matrix 
M is positive definite. Besides that of a-quartz under 
compression, there are other known examples of soft, flat 
branches. In the perovskite ferroelectrics, for example, 
there are optical branches which are unstable in the high- 
symmetry cubic structure, which,, in that case, result in 
a ferroelectric phase transition.E3 

We suppose that the subset {fj,} of the set {A} of all 
branches are soft and flat. We construct their associ- 
ated lattice Wannier functionsEj from the corresponding 
normal modes via the transformation 



Y, 



hi 



= jz E Ap X {k) exp (-ikRi)X k x 
k\e{fj,} 

= y]%,mq»ma- (10) 
ma 



Since the Sip <ma decay exponentially in \R m — Ri\, the 
Yip are local modes, optimally localized by proper choice 
of the phases and amplitudes of the unitary matrix el- 
ements Ap\{k); the flatter the branches {//}, the more 
localized. In the Einstein limit of zero dispersion, the 
Yip are confined to a single cell. Also, the Hamiltonian 
contains no harmonic coupling between Yip and Xk\> for 
those A' not contained in {fi}. These local modes form 
a basis for the construction of an effective HamiltonianE£l 
much easier to use for simulationsEZrca and, in our case, 
much easier to model than the original Hamiltonian. 



IV. THE MODEL HAMILTONIAN 

Ignoring the uniform strain e, the Hamiltonian in 
Eq. (^) can be decomposed into contributions from the 
soft flat branches {fi}, the remaining branches {//}, and 
their interaction arising from anharmonicity in V[i?]: 



H = H 



H. 



H 



WIm'}- 



(11) 



In developing the model Hamiltonian, only the contri- 
bution of the soft, flat branches, H^y, is considered ex- 
plicitly. H{^y is replaced implicitly by the Hamiltonian 
of a thermal reservoir at temperature T and #{//,}{//} is 
replaced by an appropriate coupling of the {/i} to the 
reservoir. 

In what follows we drop the subscript {^} from H^y. 
The dependence of the potential V on the local mode 
amplitudes Yip can be expressed as a cluster expansion 
so that H takes the form 



^ = £ 



p2 

2Mp 



£ ^(%) + £ Vlm{Yip,Y m p.) 



IP ' l tym 

+ V lmn(Yip,Y m p l ,Y n p,,) + 0(Y i ), (12) 

where f3, (3', /3", ... G {fi}. We now assume all essential 
physics remains if we limit the anharmonicity to the local 



term Vi(Yip), a restriction which simplifies H to 



1/3 I 

+ \ £ V h?(Yit} - Ym/sMp. - Y m p,). (13) 

If we now suppose that only a single branch becomes 
soft and flat, H further simplifies, after dropping the in- 
dex p, to 



(14) 



As described in the Sec. II, we are dealing with strongly 
first-order phase transitions. The simplest functional de- 
pendence of Vi (Yi ) which can lead to such phase transi- 
tions is that of a sixth degree polynomial, 



V, 



1 



-AY? - -BY, 



CYf, 



(15) 



where B and C are positive, and A can be positive 
or negative. Not enough is gained by including coeffi- 
cients Vim beyond nearest neighbors to warrant the added 
complexity and increase in the number of parameters. 
The harmonic coupling term is then taken to be simply 
\V\ Y^u m)0^ ~ Y m ) 2 , where the notation (l,m) implies 
that the sum is restricted to nearest-neighbors only. 

We now scale displacement by the factor L = y 'B/C, 
time by to — B 2 /MC, momentum by ML /to, and energy 
by B 3 /C 2 . The resulting scaled Hamiltonian has the 
form 

H = J2l p " + J2 + \k e («j - u m f, (i6) 

l l (l,m) 

where H is the scaled Hamiltonian and P; and Ui are 
the scaled local-mode momenta and amplitudes. V\ oc (ui) 
takes the form 



V\oc{ui) 



where 



and k in Eq. (|16|) is 



1 2 1 4 , 1 6 

—au, u, H — u, , 

2 1 4 1 6 " 



AC 



B 2 ' 



ViC_ 
B 2 ' 



(17) 



(18) 



(19) 



For dcfinitcncss we take V\, and therefore k, to be posi- 
tive. H thus depends only on the two parameters a and 
k. We take it to be classical. 

The local mode amplitudes ui evolve with time under 
the influence of H and the coupling to the reservoir. We 
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presume ui to follow Langevin dynamics with the equa- 
tion of motion 



dVtot 1. , „ 
ui = — in + Fi, 

OUl T 



(20) 



where 



Vtot = ^2 V loci u rn) + 7} k ^2 ( U ™ ~ ""^ ^ 
m ~ {rn.n) 

is the total potential. Fi is a Gaussian random force 
exerted by the reservoir with moments 



(Fi(t)) = 0, and 
(F l {t)F m (t'))=F6i m 6(t-t'). 



Equipartition requires that 
(*?> = 



-k R T 



(22) 
(23) 



(24) 



with T the reservoir temperature. The fluctuation- 
dissipation theorem requires that 



1 

T 



(25) 



Thus we are faced with a four-parameter problem 
(a, fc, t, T) even with this simplest of models. Achieving 
a detailed understanding of the physics of the model re- 
quires extensive numerical simulation. However, many of 
its important qualitative features are readily extracted, 
yielding important insights into the nature of pressure 
amorphization in those materials to which the model 
might be relevant, as described in the next four sec- 
tions. Thp-,results of the simulations will be reported 
separatelyCJ 



V. PRESSURE AMORPHIZATION 

We now proceed to eluciate the potential energy land- 
scape of our model. For a > 3/16 and Vfc in Eq.(21), 
Vtot has a global minimum at ui = 0, VZ, at which 
Vtot = 0. We assign this minimum to the low-pressure 
crystal structure C described in Sec. II. When a < 3/16 
and Vfc, there are two equivalent global minima with 



ui = w 



Uj 



1 



, V7, and 



V/ 



at which 



-N[-a 
3 L 4 



(- 



a 



,3/2 



]<0. 



(26) 
(27) 

(28) 



We take this structure to represent the high-pressure 
crystal structure C of Sec. II. Thus, there is a first-order 



phase transition at a = 3/16 from C to C as a decreases 
monotonically with pressure p irrespective of the value of 
k. The two spinodals occur at a = 1/4 and a = 0, Vfc. 
For the phase transition C «-» C to be strongly first- 
order in the sense used in Sec. II., a variation in a from 
1/4 to can be thought of as occuring via substantial 
change of pressure. The barrier between C and C at the 
transition (where a = 3/16) is 1/96, so that the reduced 
temperature T must be taken substantially lower than 
1/96 for strong hysteresis to occur, as it does in type I 
materials. 

We first consider the Einstein case in which fc = 0. 
For 3/16 < a < 1/4, 2 N degenerate local minima exist 
in Vtot with energy given by Eq.([2^), except that Vtot is 
positive in this range. The magnitude of ui is given by 
w in Eq. (|26j) , its sign being arbitrary. As a falls below 
3/16, all 2^ phases become equivalently stable relative 
to C. Among these 2 N phases, two correspond to single- 
crystal C, and an exponentially small fraction to poly- or 
micro-crystalline C"; yet the vast majority of the struc- 
tures are amorphous in that there is at most short-range 
order in the ui even though a clear topological mem- 
ory of the original crystal structure C remains. Barriers 
keep the system in C until the spinodal at a = is ap- 
proached. At the spinodal, C becomes locally unstable 
relative to all 2 N phases. The probability of the system 
finding its way from C to one of the amorphous phases is 
vastly higher than that to one of the exponentially fewer 
polycrystalline C phases. 

More generally, pressure amorphization occurs if (i) 
there is a region of the a-k plane in which the a phases 
are more stable than the C phase but metastable with 
respect to C and (ii) if r and T are such that a is pref- 
erentially accessible from C either irreversibly (types I, 
II, and IV) or temporarily (types III and V) . In the next 
three sections we shall argue that this indeed holds for 
the 2-site case, the 3-site case, and the AT-site case by ex- 
amining the topology of the potential-energy landscape 
presented by V[u]. In reading these sections, it is im- 
portant to bear in mind that fc determines how fiat the 
branch is and a how soft. A soft, flat branch has small 
a and small fc. 

Note that the previous discussion straightforwardly ap- 
plies to the case in which there are v flat branches. V[Y] 
is then defined on a i/A-dimensional configuration space. 
Around every minimum in this space there is a basin of 
attraction bounded by a separatrix, a hypersurface of di- 
mension i/N — 1. As in the single- branch case, PA will 
occur when the amorphous minima a are preferentiably 
accessible from C. 



VI. THE 2-SITE MODEL 

To begin, we study pressure amorphization in our sim- 
ple Hamiltonian through analysis of the topology and ki- 
netics of a 2-site model. In what follows the local-mode 
amplitude and momentum of the first site are x and P x , 
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FIG. 1: Contour plot of the potential V(x,y) of the 2-site 
model in the xy plane. Minima are indicated by white circles 
and saddle points by open squares; the maxima are clear by 
inspection. The larger the symbol, the lower the value of 
V{x,y). The crystalline minima are indicated by C or C' , 
the intermediate minima by i, and the amorphous minima by 
a. 



Soft mode 



c 



c 




Soft 
branch 



3/16 1/ 



FIG. 2: The a-k phase plane of the 2-site model divided into 
regions of distinct critical point structure. Maxima are in- 
dicated by open circles, minima by solid circles, and saddle 
points by open squares. The sizes of the circles reflect the rel- 
ative magnitudes of the extrema. Each region is labelled by 
the stable crystal structure within it, except for one. Within 
that region, marked PA, the critical point structure is the 
same as for k = 0, < a < 3/16, and pressure amorphization 
is possible as discussed in the text. 



those of the second site y and P y . The Hamiltonian is 
H = T + V, (29) 

where 

T=\{P 2 X +P%) (30) 

and 

V = ia(z 2 + y 2 )-i(x 4 + t /) + i(z 6 + y 6 ) 

+ \k{x-yf. (31) 

In Fig. 1, we show in the xy configuration plane the crit- 
ical points of V(x,y) for k — and < a < 3/16. C 
occurs at x = y = 0, C at x = y = ±w. There is only 
one saddle point separating the C, C", or a basins from 
each neighboring intermediate i basin. 

In Fig. 2, we divide the a-k phase plane into regions 
each of which has a distinct critical point (CP) structure, 
labelling each region (with one exception) by the stable 
crystal structure within it and displaying the CP struc- 
ture. The exception, shaded, is the region within which 
V retains at finite k for < a < 3/16 the same topol- 
ogy, i.e. CP structure, that it has for k = as shown in 
Fig. 1. It is only within this region that PA is possible 
because the a minima exist, they are stable relative to C, 



and both a and C can be reached from C only by trajec- 
tories which must pass through intermediate basins i so 
that the transition C to a can compete kinetically with 
C to C. 

For PA to occur within this 2-site model, trajectories 
from C to a must be traversed with some probability rel- 
ative to those from C to GL This happens if the following 
seven conditions are met £9 (i) The minimum must exist 
and be stable relative to C; (ii) Intermediate basins of 
attraction must exist, and all statistically significant tra- 
jectories from C to C or a must pass through them; (hi) 
For hysteresis, both V ic — V c and V ia — V a must be much 
larger than ksT at a — 3/16; (iv) For ia to be accessibile 
from ic, one requires Vi C — Vi a ^ and slow relaxation in 
the i basin; (v) For the probabilities of transitions from 
i to a and from i to C to be of comparable magnitudes, 
Via — Vic' must not be much larger than fc^T; (vi) For the 
transition from i to a to be irreversible, or likewise for 
there to be a very long resident time in a before a tran- 
sition from a to C, Vi a — V a ^ ksT and fast relaxation 
must occur in the basin; (vii) Finally, for consistency of 
(v) and (vi), t must not be too low. Detailed quali- 
tative considerations strongly suggest that these condi- 
tions can be met only for a and k in the shaded region of 
Fig. 2, provided that the remaining parameters T and r 
assume appropriate values. The actual determination of 
the probability of C — > a relative to that for C — > C for 
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various values of T and r requires numerical simulation, 
the results of which will be reported elsewhere. 

We have supposed that a is a monotonically decreas- 
ing function of pressure. For simplicity, we shall further 
suppose that k depends weakly on p as a decreases. The 
effect of pressure change on the system can be repre- 
sented by a trajectory in the a-k plane. One can readily 
see that the five types of behavior described in Sec. II 
have direct parallels in the responses of our model system 
along different trajectories in the a-k plane. For exam- 
ple, a trajectory which starts at a > 1/4 and moves well 
into the PA region before reversing and returning cor- 
responds to type I behavior, with p ca corresponding to 
< a ca < 3/16 and p ac corresponding to 3/16 < a ac < 
1/4. Type II corresponds to a trajectory that starts in 
structure C at a = «o with 3/16 < ao < 1/4, transforms 
into a below < a ac < 3/16, and can return to ao while 
still being amorphous. Trajectories terminating in the 
PA region near its upper right boundary correspond to 
behavior typical of type III materials. A trajectory which 
starts at a > 1/4 and passes through the PA region into 
that with a < would correspond to the situation in a 
type IV material. Finally, trajectories which, on their re- 
turn to larger a, pass through the upper right boundary 
can undergo an a — > C transition inside the boundary, 
remaining in C upon pressure reduction until some a in 
(3/16, 1/4) is reached; these would correspond to transi- 
tions in type V materials. 

VII. THE 3-SITE MODEL 

We now discuss an extension of the model described in 
Sec. VI to three sites. We label the third site by z and 
suppose it is coupled only to y as in a linear triatomic 
molecule. An additional term is added to T in 

Eq. (|30|) , and V becomes 

V=±a(x 2 +y 2 + z 2 )-^ + y i + z< 1 ) 

+h x * +y 6 + z e ) + h[(x - yf + {y- zf]. (32) 
o 2 

Just as in the 2-site model, there is a region in the a-k 
plane within which the CP structure and the topology 
of V remain the same as for k = 0. It is within this 
region that PA can occur. Its shape is qualitatively the 
same as that for the 2-site model, bounded by the ct-axis 
from to 3/16 and on the right by a vertical segment at 
a = 3/16. The upper boundary segments, however, are 
displaced towards smaller k. 

The k — minima are shown in Fig. 3. For small 
but finite k in the PA region, the structure is com- 
pressed perpendicular to the x = y and y = z planes. 
There are now six classes of minima: one C, six i%, 
six «2) six i' 2 , six a, and two C , using the notation of 
Fig. 3. The sequence of values of V at the minima is 
V c > V it > Vi> > V i2 > V a > Vc- The sequence at 
the saddles is V ci > V hi > 2 > V ili2 > V i2a , V^ ai V^c 




FIG. 3: Minima of the potential energy landscape of the 3-site 
model. 

and V ci - V ni > 2 > V lli ' 2 - V ni2 , etc. These sequences 
have the following implications. As before, transitions 
between the basins occur over separatrices near the tran- 
sition states (saddle points). Relaxation can be fast 
enough to bias the transformation down hill but still slow 
enough that enough energy remains to get over each suc- 
cessive barrier with little difference between rates over 
the primed and unprimed barriers. To reach C from C, 
the system must follow the sequence C — * i — s- i% — ► C 
for which there are 12 independent routes. To reach a 
from C, the system can follow either of two sequences, 
C — > i — > %2 — > a or C — > i —>■ i' 2 — > a, for which there 
are 36 routes, or 3 times as many as for C — > C . This 
increase in the number of paths from C to a relative to 
that from C to C compensates for the more favorable 
kinetics of the final step from i 2 to C , increasing the 
probability of C — > a relative to that of C — > C over 
that of the two-site model. 



VIII. THE AT-SITE MODEL 

We now jump to a model with a macroscopically large 
number of equivalent sites each of which is coupled to 
a finite number of nearest neighbors via the harmonic 
term in V. Once again, there is a domain in the a-k 
plane in which the topology of V is the same for finite k 
as for k = 0. Its lower boundary is unchanged at k = 0, 
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< a < 3/16. Its upper boundary is further compressed, 
the degree of compression increasing with the number of 
nearest neighbors. Just as the minima for the 3-site case 
lie at the center, face centers, edge centers, and corners of 
a cube, so are the minima arranged on an TV-dimensional 
hypercube in the ^-dimensional configuration space for 
k = 0. As k increases, the CP structure distorts without 
change of topology up to some k which depends on a. 
The existence of this finite domain of topological stability 
is guaranteed by the analyticity of V, its gradient, and 
its Hessian in the configuration space. There is now one 
minimum of type C, 3^ — 2^ — 1 of type i, 2 N — 2 of 
type a, and 2 of type CL, The sequence of values of 
V is V c > V% > V a > Vc'S^ The number of a minima is 
exponentially larger than the number of C minima, even 
if one extracts the polycrystalline C structures from the 
above number of a's and adds them to the C"s, as one can 
see from the following estimate. Let n be the minimum 
grain size. Then 2 N / n is a lower bound but still a good 
estimate of the number of polycrystalline structures. The 
corrected ratio of C"s to a's is thus estimated to be 



2~ JV(l-l/n) 



1 - 2- JV ( 1 - 1 /«) JV- 



-> 0, with n > 1, 



(33) 



which is still exponentially small. The number of paths 
from C to a is correspondingly exponentially larger than 
that from C to C, overcoming any kinetic advantage of 
any single path from C to C . In conclusion, the key 
to PA within this model is the existence of a region of 
the a-k plane within which the k = topology persists, 
a region of a flat (small k) and soft but not unstable 
(0 < a < 3/16) branches together with r and T in a 
suitable range. 



IX. RELATION TO CONVENTIONAL 
NUCLEATION AND GROWTH MECHANISM 

Suppose the transformation C — ► C were to occur 
by the conventional nucleation and growth mechanism 
with / the nucleation rate and R the growth rate.EJ 
For I/R finite, a polycrystal of C results. In the limit 
I/R — » 0, a single crystal results. In the limit I/R-* oo, 
a nanocrystalline material results, and we can imagine 
that in the ultimate limit in which nucleation totally 
dominates growth, an amorphous phase results when the 
crystallite is so small that there is no distinction between 
bulk and grain boundary. In the classical theory of nu- 
cleation, the critical radius of the nucleus of C is 



A 5 



(34) 



Here Ag is the bulk Gibbs free-energy difference between 
C and C, and a is the interfacial free energy. The cor- 
responding barrier to nucleation is 



AG* = 



16tt cr J 
~(A5? 



(35) 



For nucleation completely to dominate growth and for 
the critical size nucleus to be so small that the result is 
an amorphous phase, it is obvious from Eqs. ( |34] ) and fl35| ) 
that a must become small, a measures the free-energy 
cost of a rapid local structural change from C to C . 
In our simple model, the energy cost of such a structural 
change is simply proportional to the coupling constant k. 
Thus, our finding the possibility of pressure amorphiza- 
tion in the region of the a-k plane within which the k = 
topology of V is preserved is completely consistent with 
nucleation theory, the smallncss of k corresponding to 
small a and the existence of intermediate minima corre- 
sponding to finite Ag and slow growth rate. To conclude, 
PA occurs when nucleation totally dominates growth. 



X. SUMMARY AND CONCLUSIONS 

In this work, we make a clear distinction between two 
types of amorphous phases. The first, the conventional 
glassy phase, is not topologically equivalent to any crys- 
talline phase, and glass formation from the crystal is a 
reconstructive transformation^ Such glasses are formed, 
for example, by sufficiently rapid cooling from the melt. 
A conventional glass can form via pressure amorphiza- 
tion if (i) the limit of local stability of C relative to its 
melt (the spinodal) is approached, and (ii) if the glass 
transition temperature of the melt is then -higher than 
the equilibrium melting temperature of C.cru Here we 
are concerned not with conventional glasses, but with 
pressure-induced amorphous phases in which the struc- 
tural instability that leads to the amorphization is asso- 
ciated with one or more soft, flat branches of the phonon 
spectrum and in which memory of the original topology 
of the low-pressure crystal structure is retained through 
the remaining stable branches. In such cases, the amor- 
phization arises from a random displacive or orientational 
transition instead of from a reconstructive transition. 

This type of pressure amorphization is driven by an 
underlying strongly first-order structural phase transi- 
tion C — > C" under compression. In our view, as C be- 
comes thermodynamically unstable relative to C, there 
is one or more branches which flatten with increasing 
pressure. Local modes in the form of lattice Wannier 
functions provide a natural means of description of such 
flat branches. We have supposed that the nonlinearities 
associated with the first-order phase transition arise pri- 
marily in the amplitudes of the individual local modes, an 
inessential simplification. We have argued that there is 
a range of pressures within which the local nonlinearities 
and relative flatness of the relevant branch or branches 
allow an exponentially large number of amorphous struc- 
tures to be metastable and preferentially kinetically ac- 
cessible over C relative to the thermodynamically stable 
single crystal or polycrystalline C phase. Such metasta- 
bility and kinetic preference occur within a restricted re- 
gion of the parameter space defining the potential energy 
landscape in the configuration space of the local modes, 
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as shown explicitly for the 2-site model in Fig. 2. The 
ambient temperature must be low enough to lead to sub- 
stantial hysteresis (T <C 1/96), so that C persists well 
into its domain of metastability as p increases. Once 
the structural transition is initiated, it proceeds downhill 
via an exponentially large number of paths through the 
basins of attraction of intermediate structures without 
being trapped there, requiring relaxation rates in a suit- 
able range. If trapping were to occur, it would correspond 
to intermixed amorphous a, intermediate, and crystalline 
C regions, resulting in observable Bragg peaks associated 
with C. Further compression (i.e., further reduction of 
a) would result in a drop in the fraction of C, decreas- 
ing Bragg peak intensities. That is, trapping implies a 
broadened C to a transition. The intermediate basins 
correspond to a random admixture of both C-like and 
C"-like local configurations. 

Normally, first-order structural transitions are dis- 
cussed in terms of nucleation and growth. In the present 
picture, the critical size nucleus is not much more than 
a single unit cell (with minor displacements in neigh- 
boring cells) and nucleation overwhelms growth, leaving 
the system in an amorphous structure when nucleation 
is complete. 

The type of pattern observed in the structural transi- 
tions caused by pressure changes depends on the partic- 



ular trajectory through the parameter space of the po- 
tential which is induced by the changing pressure. 

The picture developed in this paper can embrace all 
of the phenomena which have been observed in associa- 
tion with pressure amorphization. Nevertheless, the dis- 
cussion has been largely qualitative. Obtaining a quan- 
titative demonstration that these simple models admit 
pressure amorphization requires numerical simulation.E2l 
Once that has been demonstrated, such simulations 
should be carried out for real materials since substan- 
tial experience has already been built up in constructing 
effective Hamiltoni ang | pf . the type considered here for 
specific substances.c3E3E3 Among the tasks for simula- 
tions is to distinguish whether the occurence of pressure 
amorphization results in a conventional glass or leaves the 
crystalline topology invariant for a particular material. 
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